Canted Antiferromagnetic Order of Imbalanced Fermi-Fermi mixtures in Optical 

Lattices by Dynamical Mean-Field Theory 
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We investigate antiferromagnetic order of repulsively interacting fermionic atoms in an optical 
lattice by means of Dynamical Mean-Field Theory (DMFT). Special attention is paid to the case of 
an imbalanced mixture. We take into account the presence of an underlying harmonic trap, both in 
a local density approximation and by performing full Real-Space DMFT calculations. We consider 
the case that the particle density in the trap center is at half filling, leading to an antiferromagnetic 
region in the center, surrounded by a Fermi liquid region at the edge. In the case of an imbalanced 
mixture, the antiferromagnetism is directed perpendicular to the ferromagnetic polarization and 
canted. We pay special attention to the boundary structure between the antiferromagnetic and the 
Fermi liquid phase. For the moderately strong interactions considered here, no Stoner instability 
toward a ferromagnetic phase is found. Phase separation is only observed for strong imbalance and 
sufficiently large repulsion. 

PACS numbers: 71.10.Fd, 75.50.Ee, 37.10. Jk, 67.85. Lm 



I. INTRODUCTION 

Ultracold atoms in optical lattices provide a versatile 
laboratory for interacting quantum many body systems^. 
One of the major challenges in this field is the experimen- 
tal investigation of quantum magnetism in atomic mix- 
tures. Impressive experimental progress in this direction 
has already been made. The first important step in ex- 
periments with fermionic atoms in optical lattices was the 
experimental observation of the Fermi surface 2 -. Recent 
experiments with bosonic atoms directly observed corre- 
lated particle tunneling^ and superexchange^, which are 
the basic mechanisms underlying quantum antiferromag- 
netism. Moreover, strong experimental evidence for the 
fermionic Mott insulator state was obtained, both by the 
local probe of observing reduced double occupancy^ and 
the global probe of observing a plateau in the cloud size 
when the system is compressed^. A recent experiment in 
a system of spin-1/2 fermions without optical lattice indi- 
cates a Stoner instability toward a ferromagnetic state for 
strong repulsion^. These are important steps on the way 
toward realization of strongly correlated many-fcrmion 
states^. Currently the experimental temperatures are 
still higher than the critical (Neel) temperature, below 
which antiferromagnetic order is predicted to develop^^. 
Most accurate theoretical estimates for the entropy per 
particle below which long-range antiferromagnetic order 
is expected yield a value of S/N « ln(2)/ 2 11 i 12 , whereas 
current experiments reach an average entropy which is 
still a factor 2 higher—. 

Ultracold atomic system offer the unique possibility 
to control the relative densities of the two spin compo- 
nents, as alrcday has been demonstrated in experiments 
without the presence of an optical lattice^^. Experi- 
mentally, the density imbalance is precisely tunable by 
means of radiofrequency sweep a 14 i 15 and stable due to 
the suppression of spin-flip scattering processes in cold 
atomic gases. This realizes an imbalanced spin mixture, 



in which the 5f/(2)-symmetry is broken by an artificial 
magnetic field. When the density of atoms corresponds 
to one particle per lattice site, the ground state of this 
system is expected to be a canted antiferromagnet, with 
antiferromagnetic order characterized by a Neel vector 
directed perpendicular to the applied field. However, 
experimentally ultracold atom systems arc always con- 
fined by an external harmonic trapping potential, which 
leads to an inhomogeneous system. If the total parti- 
cle number is sufficiently high, in the center of the trap 
a region with particle density per site close to one will 
develop, where antiferromagnetic order is stable at suf- 
ficiently low temperature o 16 i 17 . The edges of the system 
have lower filling; they are Fermi liquid regions without 
spin order. If the total particle number is even higher, 
also in the trap center a Fermi liquid with particle den- 
sity higher than one or a band insulating state can exist. 
In that case antiferromagnetic order can be stable in a 
shell around this Fermi liquid^. This poses interesting 
questions regarding the nature and the stability of spin 
order, which we will address in this paper by means of 
(Real-Space) Dynamical Mean-Field Theory. 

These issues have recently also been investigated by 
other methods. For homogeneous systems described 
by the hole-doped Hubbard model, both commensu- 
rate and incommensurate spin-density-waves have been 
predicted^— . By mapping to an effective spin model, 
the critical temperature for canted antiferromagnetic or- 
der was calculated and topological excitations of imbal- 
anced mixtures were studied 2 ^. A Hartree-Fock static 
mean-field theory for balanced mixtures in a trap pre- 
dicts that antiferromagnetism can coexist with param- 
agnetic states in various spatial patterns, for example 
antiferromagnetism in the center of the trap surrounded 
by a hole-doped atomic liquid or antiferromagnetism in a 
ring with a Fermi liquid in the center and at the edge 2 ^. 
For imbalanced mixtures, this approach predicts canted 
order perpendicular to the (artificial) magnetic field up 



to moderate values of the repulsion^ 3 -. Very recently the 
Hartree-Fock approach has also been applied to larger 
repulsion: in addition to canted antifcrromagnctism, a 
critical interaction was found, beyond which the Stoncr 
instability drives a ferromagnetic transition at the edge 
of the system, where the particle density is lower than 
half-filling^. 

A Real-Space Dynamical Mean-Field (R-DMFT) 
study of antiferromagnetism in a harmonic trap has also 
been performed, but so far without allowing for the pos- 
sibility of canted antiferromagnctic orde r 16 ' 17 . For the 
case of an imbalanced mixture, this constraint lead to 
the prediction of phase separation between the majority 
component in the center and the minority component at 
the edge for sufficiently strong repulsive interactions and 
large values of the imbalance 1 ^. 

Here we perform a full R-DMFT study, which includes 
the possibility of canted order. Unlike static Hartree- 
Fock mean-field theory, DMFT is a non-perturbative 
method which is reliable both for strong and weak in- 
teractions in sufficiently high dimensions. Local correla- 
tions are included exactly^—. R-DMFT thereby takes 
the inhomogencity induced by the presence of a harmonic 
trap into account in a fully consistent way. 

For imbalanced systems we indeed observe canted an- 
tiferromagnctic order. We consider weak to moderately 
strong interactions, for which no Stoner instability to- 
ward spontaneous ferromagnetism is found: in the case of 
a balanced mixture the wings of the systems are always 
paramagnetic. Only upon applying a finite amount of 
imbalance, the system gets polarized and ferromagnetic 
order starts to develop. We generally also do not ob- 
serve the phase-separation scenario for large imbalance. 
Instead, the canted antiferromagnetic order allows for 
a continuous transition between balanced antiferromag- 
netism order and fully imbalanced ferromagnetic order. 
Only for large values of the interaction and strong imbal- 
ance phase separation occurs. 

We compare our full R-DMFT results with calcula- 
tions based on a local density approximation in combi- 
nation with DMFT, in which the harmonic trap is in- 
corporated by a spatially varying chemical potential. As 
for the balanced cas o 16 i 17 we find that the total density 
is well approximated by the local density approximation, 
but a strong proximity effect is observed for the anti- 
ferromagnetic order: the staggered antiferromagnetism 
as obtained by the full R-DMFT calculation extends to 
regions where the local density approach predicts a para- 
magnetic solution. 



II. MODEL 

Repulsively interacting fermions in a sufficiently deep 
optical lattice are well described by the single-band Hub- 



bard Hamiltonian in the tight-binding approximation 
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where n ?(T = c icr Ci a , and Ci a (c icr ) are fermionic annihi- 
lation (creation) operators for an atom with spin a at 
site i, J is the hopping amplitude between nearest neigh- 
bor sites (ij), U > is the on-site interaction, \i a is 
the (spin-dependent) chemical potential and Vi = VqtI 
is the harmonic confinement potential. We also define 
|U = 5(jUf + fJ-i) and A/i = /i^ — /j,± as the average chemi- 
cal potential and difference in chemical potential, respec- 
tively. Although Afi acts as a magnetic field, experimen- 
tally the imbalance is not induced by a physical magnetic 
field, but by directly controlling the difference in particle 
number. The parameters of this model can be tuned in 
experiments by changing the intensity of the optical lat- 
tice and via Feshbach resonances 1 -. In the following, we 
take the lattice constant to be a = 1. 



III. METHOD 

To obtain the ground state properties of this system, 
we apply R-DMFTikiT-Sa-^ 7 .. Within R-DMFT the self- 
energy is taken to be local (which is exact in the infinite- 
dimensional limit 3 ^^) but allowed to depend on the lat- 
tice site, i.e. Y*ij a (iu) n ) = T,„'{ioj n )5ij, where Sij is a 
Kroneckcr delta. The lattice sites are described by local 
effective actions, each representing an effective Ander- 
son impurity model, which are coupled via the real-space 
Dyson equation for the Green's function. Details of the 
method have been published previously^ 6 ". 

In the present paper we use Exact Diagonalization 
(ED) 27 ' 40 ' 41 of the Anderson Hamiltonian to solve the 
local impurity actions. Within ED the spectral func- 
tion is represented by a finite number of delta peaks. 
Whereas this is sufficient for a faithful representation of 
the zero-temperature spectral function, we found it to 
lead to unphysical behavior at finite temperature, espe- 
cially away from half-filling. Therefore we restrict our- 
selves in this article to the low-temperature limit and 
only investigate ground state properties. The multigrid 
Hirsch-Fye quantum Monte Carlo method has proven 
to be a very efficient solver at finite temperatures for 
a balanced mixtur e 17 ' 37 ' 42 , but canted antiferromagnetic 
order is probably harder to obtain within this method. 
Also when using the Numerical Renormalization Group 
(NRG) methodi 6 .^ 3 -^ to solve the Anderson Hamilto- 
nian, it is problematic to describe the canted off-diagonal 
spin order. In contrast, the ED-method wc use here is 
very flexible, which also allows to incorporate off-diagonal 
canted spin order in a straightforward manner. However, 
since S z in this case is not a good quantum number for 
the individual spin components, the size of the Hilbert 
spaces to be diagonalized is significantly enlarged, which 
leads to far more time-consuming numerics compared to 
the balanced case. 
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FIG. 1: (Color online) Real-space spin profile for U/J — 10, 
A/i/J = 0.5 and Vo/J = 0.05. The vertical components 
of the arrows symbolize the local magnetization in the z- 
direction (Sf), while the horizontal components correspond 
to the (staggered) local magnetization in the x-direction {Sf). 
For all data in this paper we took the filling at the center 
equal to one, which is induced by the choice of the chemical 
potential equal to fi — U/2. To obtain the ground state ex- 
pectation values a small fictitious temperature of T/J — 0.02 
was applied. 



For very large repulsion, ED can run into unphysical 
instabilities. Therefore we consider only moderately large 
ratios of U/J here. 

In the practical implementation of ED a small but fi- 
nite temperature is used to generate the T = data, 
in order to obtain discrete Matsubara frequencies. We 
chose this temperature equal to T/J = 0.02 for the two- 
dimensional data and T/J = 0.025 for three dimensions. 



IV. RESULTS 

We now apply the R-DMFT method to spin-i 
fermions in a two-dimensional square lattice and a three- 
dimensional cubic lattice with harmonic confinement. 
We focus on the density distributions of the two species 

n la = {h la } (a =t, |), the total density n ltot = n it + m±, 



and the local spin expectation values Sf 
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in which a a (a = x,y,x) are the Pauli matrices. Here 
we have set h = 1. 

In the case of a balanced mixture the Hamiltonian is 
SC/(2)-symmetric. This means that the staggered mag- 
netization can point in any direction. We have chosen 
it to point in x-dircction. In the case of an imbalanccd 
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FIG. 2: (Color online) Radial profiles for the ground state of 
two-dimensional Fermi-Fermi mixtures for different values of 
U and A/i . Plotted are the total density (hi) = (h^ + ha/), 
and the expectation values of the spin in ^-direction (Sf ) and 
in x-direction (Sf ). The left column is for a balanced mixture 
(A/i/J = 0), the right column is an imbalanced mixture with 
Afi/J = 0.5. The other parameters are p, — U/2 and from 
top to bottom: U/J = 5, Vo/J = 0.03; U/J = 10, V /J = 
0.05; U/J = 20, Vo/J = 0.1. Here and in the following, spin 
expectation values are plotted in units of h. 



mixture, the SC/(2)-symmetry is spontaneously broken 
by the chemical potential difference, which acts as an ar- 
tificial magnetic field in the z-direction. In reaction, the 
staggered magnetization orders perpendicular to this, i.e. 
in the xy-plane. The remaining C/(l) symmetry is also in 
this case in our calculations broken by a small initial nu- 
merical perturbation, resulting in alignment of the spins 
along the x-direction, such that (S y ) = in all results 
presented here. 

In all the calculations reported here we have chosen 
p, = U/2, such that the system is at half- filling at the 
center. 



A. Results for two dimensions 

In Fig. [1] a typical real-space spin configuration is 
shown for an imbalanced two-dimensional system. Here 
we have chosen to label the spatial coordinates by x and 
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application of a finite chemical potential difference fer- 
romagnetic order is induced in the wings. Although we 
cannot establish the critical interaction above which the 
Stoner instability occurs, we thus find that this value is 
considerably shifted upwards compared with the value 
obtained within the Hartree Fock analysis, where spon- 
taneous ferromagnetism was observed for even smaller 
values of U/J than considered here^i. It is indeed well- 
known that the Hartree-Fock approximation underesti- 
mates the critical interaction for spontaneous ferromag- 
netism by more than an order of magnitude^. Within 
DMFT the dynamical screening of the local repulsion is 
fully accounted for and only for extremely large on-site 
repulsion a ferromagnetic ground state was found on the 
homogeneous cubic lattice^. However, this limit is ex- 
perimentally hard to reach, because the associated criti- 
cal temperature for spin order is very low. 



FIG. 3: (Color online) Radial distribution for the minority 
density (rui.) and the in-plane staggered order (Sf) (here plot- 
ted in absolute value) for large imbalance and different values 
of the repulsion: U/J — 10 (left column) and U/J — 20 (right 
column) and p, — U/2. Other parameters: a) Vo/J = 0.05, 
Ay,/ J = 1.2; b) Vo/J = 0.07, An/ J = 0.6; c) V /J = 0.05, 
An/ J = 1.6; d) Vo/J = 0.07, A/i/J = 0.9. 



z, such that the spin direction and spatial direction can 
be identified. Results of our R-DMFT calculation on the 
two-dimensional square lattice for the radial density and 
spin profiles are shown in Fig. [2 both for the balanced 
system and for the situation that imbalance is induced 
by a nonzero chemical potential difference A/j,. 



1. Antiferromagnetic region 

We first turn our attention to the antiferromagnetic 
region in the center, where the particle density is at half- 
filling. Our results show that imbalance reduces antifer- 
romagnetic order in the center, and tilts it out of the xy- 
plane by inducing a nonzero ferromagnetic z-component 
of the spin. The effect of imbalance becomes larger with 
increasing interactions. This can be understood from the 
fact that for large U/J the local spins in the insulating 
region interact via a Heisenberg Hamiltonian 



u = j^Y^Si-Sj-AnJ^sf. 



(2) 



<u> 



Here J ox = 2J 2 /U, such that with increasing U the ex- 
change coupling between the spins decreases and becomes 
weaker relative to the applied chemical potential differ- 
ence. 

We do not observe any sign of a Stoner instability in 
the paramagnetic wings for the values of U/J considered 
here, which would mean a spontaneous ferromagnetic po- 
larization for equal chemical potentials of the spin com- 
ponents. In contrast, we observe that only upon the 



2. Boundary structure 

We now turn our attention to the boundary structure 
between the antiferromagnetic core and the paramagnetic 
wings. The results in Fig. [2]show a local maximum of the 
ferromagnetic polarization in the z-dircction. This max- 
imum appears at the location where the antiferromag- 
netic order in the center disappears. The minority den- 
sity shows a local minimum at this point (cf. Fig. GJi,b), 
whereas the majority density has a local maximum. This 
feature was also observed in the Hartree Fock analysis^. 
The maximum appears because the antiferromagnetic or- 
der reduces the density difference of the two components 
compared to the paramagnetic situation: a smaller den- 
sity difference leads to a larger sublatticc magnetization 
and hence a lower energy. This mechanism suppresses the 
occurrence of phase separation in the trap center, which 
was found previously in the case where canted order was 
excluded^. In contrast, the possibility of canted order 
allows a continuous transition between the limiting cases 
of an antiferromagnetic phase with equal populations of 
the two species and a fully polarized ferromagnetic phase 
in which only one of the two species is still present. 

An interesting structure emerges in the paramagnetic 
outer region for large imbalance: the density of minority 
atoms shows a second maximum (cf. Fig. [3]) originat- 
ing from the strong repulsion which pushes them to the 
outside. This ring-like structure is the remnant of phase 
separation found before^. However, in this case it is 
no true phase separation, since also minority atoms are 
present in the central antiferromagnetic region. Only for 
large values of U/J in the limit of very strong imbal- 
ance, this outer ring of minority atoms survives when 
the antiferromagnetic order in the center disappears, im- 
plying true phase separation. In particular for U/J = 10 
we do not observe this scenario, in contrast to the case 
where canted antiferromagnetic order was not included^. 
When canted order is accounted for, the outer ring of 
minority atoms disappears before the antiferromagnetic 




FIG. 4: (Color online) Radial profiles for the ground state of 
three-dimensional Fermi-Fermi mixtures for U/J — 10, p, — 
U/2 and different values of Afi. Plotted are the total density 
(fii) = (h^ + hii), and the expectation values of the spin 
in z-direction {Sf ) and the absolute value of the expectation 
value in ^-direction \{Sf )|. The points denote results of the 
full R-DMFT calculation, whereas the solid lines are obtained 
within the LDA (TFA) approximation combined with DMFT. 
The values for A/j are: a) Aji/J = 0; b) Aji/J — 0.4; c) 
Afi/J = 0.8; d) Afi/J = 1.2; e) Ay,/ J = 1.6; f) Ay/ J = 2. 
The trap parameter is chosen as Vo/J = 0.15. 



order in the center vanishes, as visible in the data in Fig. 
|3£i) and c). For U/J — 20 still phase separation occurs, 
as shown Fig. [3]b) and d): for large imbalance the anti- 
ferromagnetic order in the center breaks down, but the 
ring of minority atoms surrounding the phase separated 
central region with only majority atoms is still present. 
Note that this phase separation cannot be identified with 
the Stoncr instability, because it only happens for large 
chemical potential difference. Moverover, it means that 
the complete central region, including the part where the 
density is at half-filling, is fully polarized, whereas the mi- 
nority atoms are located in a shell around this at low ma- 
jority density. In contrast, the Stoner instability favors a 
scenario, where the region at half-filling still supports an- 
tiferromagnetic order, whereas the edges are completely 
polarized^!. 



B. Results for three dimensions 

Results for the ground state (T/J — > 0) in three di- 
mensions are presented in Fig. 2] We compare the full R- 
DMFT results with a local density approximation (LDA) 
(previously also denoted as a Thomas-Fermi approxima- 
tion (TFA)i£), where the trap is modeled within DMFT 
by a locally varying chemical potential combined with the 
homogeneous density of states of the cubic lattice. In or- 
der to facilitate the comparison, the absolute value of the 
staggered spin expectation value in x-dircction is plotted 
in Fig. SI This makes it impossible to distinguish ferro- 
magnetic and antifcrromagnctic order in the figure; how- 
ever, as in Fig. [2]thc in-plane xy-spin order described by 
S x is always staggered. For the total density we observe 
good agreement between the LDA+DMFT and the full 
R-DMFT results; the only difference being that the LDA 
results show a small discontinuity, which is smoothened 
in the R-DMFT calculation. In contrast, the agreement 
between the LDA and R-DMFT results for the spin or- 
der parameters is far less good, as also observed for the 
case of balanced mixture a 16 ' 17 . In particular, antifcrro- 
magnetic spin order extends much further into the region 
with total density lower than half-filling (i.e. rii < 1) 
than LDA analysis predicts. This is a consequence of a 
proximity effect of the antifcrromagnctic insulator close 
to the paramagnetic boundary layer. 

Examining the numerical results in more detail we ob- 
serve similar features as for the two-dimensional case in 
Fig. [2] A maximum of the ferromagnetic polarization 
in the z-direction is also visible in the data in Fig. 21 
although it is more pronounced in the LDA curve than 
in the R-DMFT data, and in general less pronounced 
than in two dimensions. This is because the ratio of the 
interaction to the band-width U/2zJ (z being the num- 
ber of neighbors) is chosen smaller here. This is also 
the reason that phase separation is not observed for the 
case of strong imbalance, even though large values of A/i 
are considered. The minority atoms still form a shell for 
large imbalance, but this shell disappears before the an- 
tifcrromagnctic order in the center is destroyed. In order 
to observe phase separation, a stronger repulsion would 
be needed. Also, a Stoner instability is not observed in 
the three dimensional case for the value of the repulsion 
chosen in this case. 

Our data in three dimensions arc summarized in Fig. 
\E\ Here we plot global obscrvablcs as a function of 
the chemical potential difference: the total magnetiza- 
tion in x and z direction S^ = ^ZjWf'*) an d the 
staggered magnetization in x and z direction iS^L = 

Si( — l) lx+ly (<§f' z ), normalized to the total particle num- 
ber iVtot = 2i(^it" + foil)- Fig [S] shows that the total 
magnetization always points in the z-direction; which is 
also the direction in which it is induced experimentally 
via the spin imbalance. This means that the system is 
not frustrated and, for the parameters considered here, 
the interesting scenarios like the one pursued in Ref.— 
are not realized. 




FIG. 5: (Color online) Global observables for the three- 
dimensional system obtained by R-DMFT at U/J = 10, 
jl — [7/2 and Vo/J = 0.15 as a function of the chemical 
potential difference A/i. Shown are the total magnetization 
in x and ^-direction S^ and the staggered magnetization in 
x and z-direction S^j.'* ' , all normalized to the particle number 

/Vtot 



CONCLUSIONS 



interesting boundary structure to the paramagnetic edge, 
which we investigated in detail. 

Due to the possibility of canting, the antiferromag- 
nctism turns out to be very stable against imbalance: 
only for sufficiently large values of the repulsion and at 
large imbalance phase separation occurs. In this regime 
the central region becomes completely polarized: only 
the majority component is present and forms a band in- 
sulator in the center. The minority atoms organize them- 
selves in a shell around this central plateau. 

We do not observe a Stoner instability toward sponta- 
neous ferromagnetic order for the moderately large values 
of the repulsion considered here. The critical interaction 
to observe this phenomenon is thus necessarily relatively 
large, leading to a small critical temperature for ferro- 
magnetic order. 

Canted antiferromagnetic order is more challenging 
to detect experimentally than antiferromagnetic order 
in the z-direction, because many proposed detection 
schemes arc only sensitive to staggered order which is 
diagonal in the basis of the physical particles constitut- 
ing the system. This limitation does, however, not apply 
to Bragg scattering^, which is therefore the experimen- 
tal method of choice to detect canted antiferromagnetic 
order. 



We have solved the fermionic Hubbard model with re- 
pulsive interactions in a harmonic trap by means of Real- 
Space Dynamical Mean-Field Theory, with a focus on the 
ground state properties of imbalanced spin populations. 
We considered systems with a density of one particle per 
lattice site in the center, for which a shell structure ap- 
pears: an insulating regime in the center, in which anti- 
ferromagnetic order can arise, is surrounded by a Fermi 
liquid regime without magnetic order. Imbalance leads 
to canted antiferromagnetic order in the center, with an 



Acknowledgments 

We thank I. Bloch, E. Dcmler, W. Ketterle, A. Koetsier 
and H.T.C. Stoof for useful discussions. This work was 
supported by the Ncdcrlandsc Organisatie voor Weten- 
schappclijk Onderzoek (NWO) and the German Science 
Foundation DFG via Forschcrgruppc FOR 801, Sondcr- 
forschungsbereich SFB-TRR 49 and the DIP project BL 
574/10-1. 



I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

M. Kohl, H. Moritz, T. Stoferle, K. Giinter, and T. 
Esslinger, Phys. Rev. Lett. 94, 080403 (2005). 
S. Foiling, S. Trotzky, P. Cheinet, M. Feld, R. Saers, A. 
Widera, T. Miiller, and I. Bloch, Nature 448, 1029 (2007). 
S. Trotzky, P. Cheinet P, S. Foiling, M. Feld, U. Schnor- 
rberger, A.M. Rey, A. Polkovnikov, E.A. Demler, M.D. 
Lukin, and I. Bloch, Science 319, 295 (2008). 
R. Jordens, N. Strohmaier, K. Giinter, H. Moritz, and T. 
Esslinger, Nature 455, 204 (2008). 

U. Schneider, L. Hackermiiller, S. Will, Th. Best, I. Bloch, 
T. A. Costi, R. W. Helmes, D. Rasch, and A. Rosch, Sci- 
ence 322, 1520 (2008). 

G.-B. Jo, Y.-R. Lee, J.-H. Choi, C. A. Christensen, T. H. 
Kim, J. H. Thywissen, D. E. Pritchard, and W. Ketterle, 
Science 325, 1521 (2009). 

W. Hofstetter, J.I. Cirac, P. Zoller, E. Demler and M.D. 
Lukin, Phys. Rev. Lett. 89, 220407 (2002). 



F. Werner, O. Parcollet, A. Georges, and S.R. Hassan, 
Phys. Rev. Lett. 95, 056401 (2005). 

A. Koetsier, RA. Duine, I. Bloch, and H.T.C. Stoof, Phys. 
Rev. A 77, 023623 (2008). 

L. De Leo, C. Kollath, A. Georges, M. Ferrero, and O. 
Parcollet, Phys. Rev. Lett. 101, 210403 (2008). 
S. Wessel, Phys. Rev. B 81, 052405 (2010). 
R. Jordens et al., Phys. Rev. Lett. 104, 180401 (2010). 
M.W. Zwierlein, A. Schirotzek, C.H. Schunck, and W, Ket- 
terle, Science 311, 492 (2006). 

G. B. Partridge, W. Li, R. I. Kamar, Y. Liao, and R. G. 
Hulet, Science 311, 503 (2006). 

M. Snoek, I. Titvinidze, C. T6ke, K. Byczuk, and W. Hof- 
stetter, New J. Phys. 10, 093008 (2008). 
E.V. Gorelik, I. Titvinidze, W. Hofstetter, M. Snoek and 
N. Bhimer, Phys. Rev. Lett. 105, 065301 (2010). 
B.I. Shraiman and ED. Siggia, Phys. Rev. Lett. 62, 1564 
(1989). 
H.J. Schulz, Phys. Rev. Lett. 64, 1445 (1990). 



20 J.K. Freericks and M. Jarrell, Phys. Rev. Lett. 74, 186 
(1995). 

21 A. Koetsier, F. van Liere, and H. T. C. Stoof, Phys. Rev. 
A 81, 023628 (2010). 

22 B.M. Andersen and G.M. Bruun, Phys. Rev. A 76, 041602 
(2007). 

23 T. Gottwald and P.G.J, van Dongen, Phys. Rev. A 80, 
033603 (2009). 

24 B. Wunsch, L. Fritz, N.T. Zinner, E. Manousakis, and E. 
Demler, Phys. Rev. A 81, 013616 (2010). 

25 D. Vollhardt Correlated Electron Systems, vol 9, ed. V. J. 
Emery, (Singapore: World Scientific) p 57 (1993). 

26 Th. Pruschke, M. Jarrell, and J.K. Freericks, Adv. in Phys. 
44, 187 (1995). 

27 A. Georges, G. Kotliar, W. Krauth, and M.J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

28 G. Kotliar and D. Vollhardt, Phys. Today 57(3), 53 (2004). 

29 V. Dobrosavljevic and G. Kotliar, Phil. Trans. R. Soc. 
Lond. A 356, 57 (1998). 

30 M. Potthoff and W. Nolting, Phys. Rev. B 59, 2549 (1999). 

31 M.-T. Tran, Phys. Rev. B 73, 205110 (2006). 

32 M.-T. Tran, Phys. Rev. B 76, 245122 (2007). 

33 Y. Song, R. Wortis and W. A. Atkinson, Phys. Rev. B 77, 
054202 (2008). 

34 R.W. Helmes, T.A. Costi, and A. Rosch, Phys. Rev. Lett. 
100, 056403 (2008). 

35 A. Koga, T. Higashiyama, K. Inaba, S. Suga, and N. 



Kawakami, Phys. Rev. A 79, 013607 (2009). 

36 D.-H. Kim, J. J. Kinmmen, J. -P. Martikainen, and P. 
Torma, e-print larXiv:1009.5676l 

37 N. Bliimer and E. V. Gorelik, Comput. Phys. Commun. 
118, 115 (2011). 

38 W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 
(1989). 

39 E. Muller-Hartmann, Z. Phys. B 76, 211 (1989). 

40 M. Caffarel and W. Krauth, Phys. Rev. Lett. 72, 1545 
(1994). 

41 Q. Si, M. J. Rozenberg, G. Kotliar, and A. E. Ruckenstein, 
Phys. Rev. Lett. 72, 2761 (1994). 

42 E. V. Gorelik and N. Bliimer, Phys. Rev. A 80, 051602(R) 
(2009). 

43 G. Wilson, Rev. Mod. Phys. 47, 773 (1975). 

44 R. Bulla, T. Costi, and Th. Pruschke, Rev. Mod. Phys. 80, 
395 (2008). 

45 D. Vollhardt, N. Bliimer, K. Held, and M. Kollar, in Pro- 
ceedings of the Workshop on Ground-State and Finite- 
Temperature Band Ferromagnetism, Berlin /Wandlitz, 4~ 
6 October 2000, Lecture Notes in Physics, Vol. 580, 
(Springer, Berlin, 2001), pp. 191-207. 

46 Th. Obermeier, Th. Pruschke, and J. Keller, Phys. Rev. B 
56, R8479 (1997). 

47 T. A. Corcovilos, S. K. Baur, J. M. Hitchcock, E. J. 
Mueller, and R. G. Hulet, Phys. Rev. A 81, 013415 (2010). 



